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Abstract 



We study the role of the brane-induced graviton kinetic term in theories 
with large extra dimensions. In five dimensions we construct a model with 
a TeV-scale fundamental Planck mass and a flat extra dimension the size of 
which can be astronomically large. 4D gravity on the brane is mediated by a 
massless zero-mode, whereas the couplings of the heavy Kaluza-Klein modes 
to ordinary matter are suppressed. The model can manifest itself through the 
predicted deviations from Einstein theory in long distance precision measure- 
ments of the planetary orbits. The bulk states can be a rather exotic form of 
dark matter, which at sub-solar distances interact via strong 5D gravitational 
force. We show that the induced term changes dramatically the phenomenol- 
ogy of sub-millimeter extra dimensions. For instance, high-energy constraints 
from star cooling or cosmology can be substantially relaxed. 
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I. INTRODUCTION 



It has been known for some time that in theories with large extra dimensions the 
scale at which gravity becomes strong, M, can be much lower than the four- dimensional 
Planck mass Mp. In this approach the standard hierarchy problem takes a different meaning 
since the field theory cut-off is lowered. The four-dimensional gravity is weak due to the fact 
that gravitational flux spreads into extra dimensions. The relation between fundamental 
and observed Planck scales takes the form: 

Ml = M2+^ i?^ , (1.1) 

where R is the circumference of the extra dimensions. In this framework the Standard 
Model particles are localized on a brane whereas gravity propagates in the bulk of N new 
dimensions. Until now, it was assumed that if M ~ few TeV the = 1 case is ruled 
out by observations since it would require R ~ 10^^ cm (the solar system size). As to 
the higher codimensions, they allow R to be around a millimeter or smaller. The purpose 
of the present paper is to argue that the framework of compact flat extra dimensions can 
be dramatically modified due to the presence of the brane-induced kinetic term for higher 
dimensional graviton [^,^ 

^ind OC J d^x \J\deig\ R{x) , (1.2) 

where g is the higher dimensional metric evaluated at the position of the brane, and R{x) 
is the corresponding four dimensional Ricci scalar. As a result of this term the fiat extra 
dimension can be astronomically large, and many high energy constraints can be lifted. 

In fact, in the approach of [0,0, the extra dimension was taken to be truly infinite. In 
spite of this, it was shown that the model reproduces four-dimensional gravity due to the 
four-dimensional curvature term in the brane worldvolume theory. 

To be more precise let us consider 5-dimensional Minkowski space with a standard bulk 
gravitational action 

5bulk = / d'+'X ^\C{Gab, nABCD, $) , (1.3) 

where the capital Latin indices run over D = {A + l)-dimensional space-time. Gab denotes 
the metric of 5-dimensional space-time, TZ abcd is the 5-dimensional Riemann tensor and $ 
collectively denotes other fields. Suppose that there is a 3-brane in this space. The 3-brane 
can be realized as a soliton of the corresponding field equations. We split the coordinates in 
5-dimensions as follows: 

= {x^, y) , (1.4) 

where Greek indices run over four- dimensional brane worldvolume, fi = 0, 1,2,3 , and y is 
the coordinate transverse to the brane. In order to reduce our discussion to its main point 
the brane will be taken to have zero widthQ. In this approximation the brane action takes 
the form: 



^This is a good approximation when the mass of the field (presumably a scalar) out of which the 
brane is made is bigger than the energy scale of the low-energy effective theory. 
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•S's-brane = -T d^x \/\detg\ , (1.5) 



where T stands for the brane tension and (jf^i, = d^X^dyX^GAB denotes the induced metric 
on a brane0. In most of this work, unless otherwise stated, we neglect brane fluctuations. 
Therefore, the induced metric can be written as follows: 

9t,v {x) = G^y{x,y = 0) . (1.6) 

In general, there could be localized matter fields on the brane worldvolume. These can be 
taken into account by writing the following action for the brane: 

S^wle = ^3-brane + / d'x ,/\d^\ C{<j>) , (1.7) 

where (j) collectively denotes all the localized fields for which the four- dimensional Lagrangian 
density is C. In the classical theory, which we are discussing so far, the 4D Ricci scalar on 
the brane worldvolume is not present. Thus, the localized particles separated at a distance 
r on a brane interact via the (4 + l)-dimensional gravitational force-law, that is F ~ l/r^~^^. 
This holds as long as the classical theory is concerned. However, in the full quantum theory 
the 4D Ricci scalar will be generated (along with other terms) on the brane worldvolume. 
This is due to quantum loops of the matter fields which are localized on the brane ■ As 
a result, the following worldvolume terms should be included when one considers the full 
quantum theory: 

Sind = J V^d^ [A + R{x) + O {r^)] , (1.8) 

where M = is some parameter which depends on the details of the worldvolume 

model 0,0. A in ( |1.8| ) is an induced four-dimensional cosmological constant. The role of 
this term is to renormalize the brane tension. Furthermore, R{x) is the four- dimensional 
Ricci scalar which is constructed out of the induced metric g^v{x) defined in ( |1.6| ). In 
five-dimensional Minkowski space a brane with nonzero tension inflates Therefore, to 

avoid the worldvolume inflation we fine tune the brane tension T and the brane worldvolume 
cosmological constant A so that the net tension is vanishing 

T' = T-AM^ = . (1.9) 

This is a usual fine tuning of the four- dimensional cosmological constant. 

The graviton propagator resulting from such a system is quite peculiar. Ignoring the 
tensor structure for a moment we obtain for the corresponding Green's function the following 
expression: 

GR{p,y) = —— — ^ _2 ^ exp{-p\y\} . (1.10) 
2M3p + M p2 



^Below we will imply that the Gibbons-Hawking surface term is also present on the worldvolume 
so that it provides the correct bulk Einstein equations. 
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Here is a four-dimensional Euclidean momentum and p = For sources which are 

localized on the brane {y = 0), this propagator reduces to a massless four- dimensional 
Green's function 

GR{p,y = 0) oc 4 ' (1-11) 

provided that p ^ 1/tc- Thus, at distances r ^ rc we observe the correct Newtonian 
behavior of the potential 

V{r) oc . (1.12) 

At large distances, r 3> Tc, however, the behavior of the Green's function changes 

GRip,y = 0) oc - . (1.13) 
p 

This gives rise to a Newtonian potential which scales in accordance with the laws of a 
five-dimensional theory | 

oc ^ . (1.14) 

This somewhat puzzling behavior can be understood in two equivalent ways which we briefly 
discuss. First let us adopt the five-dimensional point of view. In this language, although 
there is no localized massless particle, there exists a localized resonance state in the spec- 
trum. The lifetime of this resonance is ~ Vc- The resonance decays into the continuum 
of modes. This can be manifestly seen using the Kallen-Lehmann representation for the 
Green's function 

Gn{p,y = 0) = ,,.3 ^\.3 , = rM^, (1.15) 
2M-^p + M-^TcP Jo s + p'^ 

where the spectral density as a function of the Mandelstam variable s takes the form: 

Pis) oc ^ -^i^ . (1.16) 

a/s 4 + s 

As Tc — >■ oo the spectral density tends to the Dirac function p{s) — >■ const.6{s) describing a 
stable massless graviton (this corresponds to the limit when the bulk kinetic term can be 
neglected). Summarizing, at distances r < the resonance mimics the massless exchange, 
and therefore mediates the 1/r^ force. At larger distances, however, it decays into the 
continuum, and as a result the force-law becomes that of a five-dimensional theory, ~ 1/r^. 



^ Transition to five dimensional regime at long distances is also a characteristic feature of the 
models with "quasi-localized" gravity [0-|l^]. However, these models usually include negative norm 
states |-|ll[. 
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A different but equivalent way to understand the above result is to adopt the point of 
view of the four-dimensional mode expansion. The analysis of the linearized equation for 
the small fluctuations shows (see Appendix) that there is a continuum of 4D massive states 
with wave-function profiles that, at the location of the brane, are suppressed by the following 
factor 

l0™(y = O)P = . J , , ■ (1.17) 

The Newtonian force on the brane is mediated by the exchange of all these Kaluza-Klein 
(KK) modes. These give rise to the expression: 

1 dm e~™' 

At any distance r the dominant contribution comes from the modes lighter than m = 1/r. 
The modes with m < l/r^ have unsuppressed wave-functions on the brane. Therefore, for 
r > Tc the interaction picture is similar to that of a five-dimensional theory. In contrast, 
when r < Vc the picture changes since the modes with m > l/r^ have suppressed couplings. 
Although the number of the modes which participate in the exchange at a given distance 
r < Tc is the same as in the five-dimensional picture, their contributions are suppressed. 
Thus, the number of the light modes effectively contributing to the exchange "freezes-out" 
and the resulting behavior of the potential is 1/r. 

As pointed out in 0, this model exhibits the van Dam-Veltman-Zakharov (vDVZ) dis- 



continuity 0,11^] in the tensor structure of the graviton propagator in the lowest tree-level 
approximation (as would any fiat space ghost-free theory where AD gravity is mimicked by 
the exchange of (a continuum of) massive spin-2 particles 0). Another disadvantage of this 
model is the low value of the crossover scale — 10^^ cm. This would imply that Newtonian 
gravity is modified at this distance, which is certainly unacceptable. 

In the present paper we show that both aforementioned difficulties are circumvented if 
the fifth dimension is compactified on a circle of macroscopic size R. We will argue that R 
can be of astronomical size. In this framework, the four- dimensional Planck scale emerges 
as 

M|. = { R + r,) . (1.19) 

Gravity on the brane is dominated by a lowest lying massless KK zero mode, which has an 
unsuppressed wave-function. As a result there is no vDVZ discontinuity problem. Contribu- 
tions from the massive KK modes give rise to small corrections to Newtonian gravity. These 
corrections are mostly constrained by the data on the motion of the inner planets. Moon 
and satellite experiments. We study constraints coming both from astronomy (planetary 
motion) and from astrophysics (star cooling), and show that the model is consistent with 
all the observations. 

II. THE MODEL 

We would like to consider a simplest possible model with a single brane embedded in 
five-dimensional space with one compact dimension. The action of our model takes the form: 
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s 



R 



dy VG 7^(5) + M 



d'^x 



R 



(2.1) 



As we mentioned before, the four-dimensional Planck mass is set by the sum of the two 
distance scales and R through the relation ( |1.19|) . For the solution of the hierarchy 
problem along the lines of Ref. |1| we can put M < TeV. This sets the largest of the two 
scales, R and rc, to be larger then ~ lO^^cm or so. It is tempting to assume that the origins 
of these two scales are related. Some possible ideas along this direction will be discussed 
later. For the time being we shall take them as the parameters of the theory. 

We are interested in the Newtonian gravitational interaction of the brane-localized 
sources. In this section we will ignore the tensor structure of the graviton propagator and 
concentrate on its scalar part. A straightforward calculation (see Appendix) leads to the 
following result for the Euclidean two-point Green's function: 



A(p,0) = - 



pi 



1 



1 + (l/pr,)2tanh(pi?/2) 



(2.2) 



For Tf.^ R this propagator exhibits the l/p^- behavior both for large and for small momenta. 



A(p,0) 
A(p,0) 



1 

p2 



1 - 
1 - 



pr 
R- 



pR^ 1 (short distances), 
pR (large distances). 



(2.3) 



Therefore the potential is four-dimensional both for large and for small distances. Note that 
the maximal deviation in the coefficient of 1/p^ is of order R/vc- 

Let us now turn to Minkowski space-time by rotating p — * ip^- The propagator takes 
the form: 



1 



A(Pm,0) = 



1 



l + {l/pMrc)2tan{pMR/2) 



(2.4) 



This expression has poles at r^PM = — 2 tan(pA,/-R/2). In order to understand better this 
result let us turn to the KK expansion. The mode decomposition is given in Appendix. 
As one would expect, there is a discrete tower of KK states. For Tc -C i? the KK masses 
approach the usual KK spectrum, m„ = 27m/ R. For ^ R all modes (except the zero 
mode) tend to the asymptotes of the tangent function: mo = 0,m„ ~ (2n — Vjn/R. This 



is in agreement with the pole structure of the propagator ( |2.4] ). Note that the level spacing 
for fixed R does not change. For arbitrary Tc the mass of the n'th state is in the interval 
((2ra - l)7r/i?, 2n'K/R). 

Let us calculate the static potential between two objects on the brane which are separated 
by a distance r. For this we perform the KK decomposition of the corresponding 5D fields. 
This is done in Appendix. The result can be summarized as follows: 



V{r) 



Go 



(2.5) 



n=l 



where for the case ^ R we can use 
find (see Appendix): 



{2n — l)7c/R and after some simplifications we 
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Gn 



R 



r,M3 \rj [{2n - l)7r/2]2 + R/r^ + {R/r^Y 
2\ /R 1 



(2.7) 



Finally, we derive the potential: 



- ^ 1 + • (2-8) 



The term in the parenthesis in this expression changes slowly with r and the magnitude of 
the change is of the order R/tc- One could think of ( |2.8| ) as being the Newton potential with 
an effective space-dependent Newton's constant. To estimate the variation of the effective 
gravitational constant we evaluate the difference of the potentials at two distinct points. 
For a given distance r in the KK sum we can neglect the modes with n > R/r (since these 
are exponentially suppressed) and replace the exponential by 1 in the remaining terms for 
which n < R/r. The variation of the Newton constant for two points ri and r2 (ri > r2) 
evaluates as follows: 

G{n)-G{r,) R 1 

G 1. -2- (2-9) 

' c n=R/ri 

Moreover, the potential can be evaluated exactly and the expression can be expanded for 
small and large values of r/R. For r/i? -C 1 and R <^ r^, we find: 



V{r) ~ ^ 



^ +-fl + -^(ln;, + C)+0 ((i?/r,)2, (r/Rf 



1 + R/vc TcK IT R\ R 



(2.10) 



where G = ln(7r/2e) ~ —0.54. Likewise, for r/i? ^ 1 and R -C Tc, the potential reads as 
follows: 



+ (^- + 0((i?/rJ^))e--/« + ... 



(2.11) 



1 + R/tc Vvr^ Tc 

The overall relative change in Newton's constant for r in the interval [0, oo] is the ratio R/r^. 



III. CONSTRAINTS 

A. Tensor Structure of the Propagator 

In the present model 4D-gravity on the brane is mediated by the zero mode graviton 
with an admixture of the tower of massive KK states. The latter contribution is suppressed 
both by the parameter a = R/vc and by the KK number n^. In order to determine the 
tensor structure of the graviton propagator let us compare how it changes when one turns 
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from a massless graviton to massive one. In 4D space the momentum- independent part of 
the tensor structure takes the form: 

l(^/^"^^/3 + _ IryM^r^"'^ . (3.1) 

This can be seen by using, for instance, the harmonic gauge for the fluctuations of the 
massless gravitational field: 

a-e^ = \ hf' . (3.2) 

On the other hand, for a massive graviton one finds: 

\[jf'-rf^ + r^^^r^^") - \^^rff^ . (3.3) 
2 3 

In the latter case, the equations of motion for a massive graviton (we use the Pauli-Fierz 
mass term) give rise to the following relation for the fluctuations: 

d^h^;::) = hf^^ . (3.4) 

Using this condition in the Einstein equations naturally leads to the tensor structure given 
in ( |3.3| )p|. Let us now see what happens in the 5D model. Before compactification one can 
choose the harmonic gauge in the bulk space: 

d'^hAB = \dBh%. (3.5) 



This leads to the tensor structure of the form ( |3.3|) . Let us consider what happens upon 
compactification. For the massive KK modes, m 7^ 0, the {^v} components of the gauge 
condition ( p.5| ) will turn into the following expression: 

d'h^;:^ = \{d^ hf^' + 9. ht^') , (3.6) 

while the {55} component takes the form: 

h^Hx) = ht^\x) . (3.7) 

Substituting ( p.7|) into (|3.6| ) we derive (|3^ ). Using this, we obtain the following expression 
for the tensor structure of massive KK gravitons: 

h^;J{p,0)f^'' ^ {f^,f'^''-^f:f'f^, (3.8) 



^ Note that the condition (3.4) would be an inconsistent gauge choice for the massless gravitons. 
This is in contrast with the vector field case, where the equation of motion for a massive vector 
field (Proca field) would require the relation d^A^^ = and this latter is also an acceptable 
gauge choice for a massless gauge field. 
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where T^i^, T'^^ are Fourier transformed matter source energy-momentum tensors. Therefore 
in the approximation R/r^ -C 1 the expression for the propagator takes the form (the tensor 
structure of the zero-mode is that of a 4D theory): 



where 



\ 2 2 / J) 

2/?°°l/l 1 \ 1 



+ 



TT^ Tc ^2 V2 



(3.9) 



(3.10) 



Thus, the lowest massive KK modes are suppressed by the ratio R/r^ = a. These contribu- 
tions might be important. Indeed, as shown in any ghost-free theory in which massive 
spin-2 KK contribution to 4D gravity is essential should suffer from vDVZ-discontinuity 
problem The reason is that in flat space any massive spin-2 states have 5 polar- 



izations. Three of these polarizations couple to the conserved stress energy tensor with the 
same strength and contribute to the gravitational potential. This induces a finite deviation 
from the predictions of Einstein gravity for arbitrarily small mass Therefore, theo- 



ries with massive graviton exchange encounter phenomenological difficulties and are severely 
constraint. We shall try to quantify the implications of this constraint in our model. 

As we discussed above, the massive modes give rise to the distinctive tensor structure 
in the propagator which alters the predictions for light bending by Sun and the precession 
of Mercury perihelion. Therefore, the massive KK modes must be adequately suppressed. 
The sum in (|3.9|) is rapidly convergent. As a result, the dominant contribution comes from 
a few lightest KK modes. Using the present experimental data , this gives the following 
constraint: 



R 

a = — < 3 X 10" 



(3.11) 



B. Planetary Motion 

We shall now examine the constraints on the size of the extra dimension R, coming from 
solar system dynamics. An object orbiting in the gravitational field will exchange a certain 
number of the massive KK modes. This is determined by the distance of the object from 
the center of the potential. For example, Pluto will effectively exchange less KK modes than 
Mercury. This will be seen as an effective change of Newton's constant and will constraint 
the possible values of R. In the present model, the total change in Newton's constant felt by 
an object falling from infinity towards the center of the potential would be R/vc- From the 
expansions of the potential (|2.10| ) and (|2.11|) it is obvious that most of the change happens 



in the range < r < R. 

Experimental bounds on the spatial variation of Newton's constant over planetary dis- 
tances are quite strong. They come from precise tracking of planetary orbits with the help 
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of radar beams and space probes. The results of those measurements are presented in the 
work [0 (an updated and detailed discussion can be found in |]18|). The modification of 
the Newton's potential assumed in is 

^V{r) = aCj-^, (3.12) 
r 

where a and A are parameters that describe deviations from standard potential. A nonzero 
value of these parameters would result in deviations from the Kepler's third law and would 
induce anomalous precessions of the planets' perihelia. The precise measurements of orbits 
of the inner planets and satellites put bounds on the values of a and A that are summarized 
in Ref . ITTlJii 



Let us relate approximately the parameters of our model with the parameters in ( p. 121 ) 
For simplicity we will use the following estimate 



Ay(r 




2 R °° g-(2n-l)7rr/i? 

(3.13) 



With this estimate we can relate the parameters of our model with the parameters in (|3.12|) 



R 

a , A~i?. (3.14) 



The bound on a is very stringent {a < 10^^, 0,|ll]) for the values of A in the range 
from the Earth radius to the Solar System size. If we assume M ~TeV we conclude that 
the size of extra dimension R must be either smaller than the Earth radius, or bigger than 
the solar system size. 

Let us discuss these two possibilities separately. We start with the case when R is larger 
then the solar system size. In this case, the estimate ( |3.14| ) is actually not the most precise 



one. The bounds on a and A shown in |17] were derived using the spatial variation of 



Newton's constant that appears in the force, rather than in the potential. We can compare 
the correction to the force in our model and the same correction for the pure Yukawa type 
modification , both shown in Fig.l . We see that the force correction in our model is a much 
steeper function on small values of r/R. 

The measured quantity is the parameter r]p that describes deviations from third Kepler's 
law. The ratio of semi-major axis can be measured (a^*^"*) and compared to prediction of 
Kepler's law (op) . The following equation relates 1] and the value of Newton's constant at 
the Earth {Ge) and at the planet's orbit {Gp). 

ap ( GpM, \"^ 



Newton's constant is the one that appears in the force law, and is the mass of the Sun. 
The bound on the parameters rjp is of the order 10^^ — 10^^° for the inner planets. 
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1 2 3 4 5 

Figure 1: The correction to Newton's force for our model (solid line, R/rc — 10^^) and for the Yukawa 
type correction (dashed line, a — 10~^). The relative correction (normalized to the strength of Newton's 
force) is given is units of 10^''. The x axis measures r/R for our model and r/A for the Yukawa correction. 



The most precise bound on the spatial variation of Newton's constant {rjp) comes from 
the relation between the Earth and inner planets, i.e., in the regime of small r/R, for R 
greater then the solar system size. In that regime, the change in force will be bigger for 
our model than for the Yukawa type correction and we have to compensate that effect by 
choosing somewhat larger values of R and r^. Numerical analysis shows that with the choice 
of i? ~ lO^^m, and R/tc ~ 10"'^, the bounds can be safely satisfied. Since the Planck scale 
is roughly given by rcM^, an increase of Tc by some factor would lower the fundamental 
Planck scale for a much milder value of a cubic root. With the new bound oi R > 10 ^^m 
and R/tc > 10^^ the fundamental Planck scale is pushed to approximately 100 GeV. 

Similarly, we calculated the contribution to the anomalous precession of Mercury and 
Mars in our model, and we found that it can be within the accepted values if the above 
limits on R and R/tc are satisfied 

The second possibility is to take the size of the extra dimension to be anywhere less than 
the Earth radius (~ 10^ m), in order to satisfy constraints from Earth- Moon and artificial 
satellite experiments. This "regime" is much safer and the bounds are very weak, as can 
be seen in [0,|l5l- In the case R ~ 10'^ m, can range form ~ 10^^ m to ~ 10^^ m, 
corresponding to a fundamental Plank scale of ~ 10 TeV down to 100 GeV, respectively. Of 
course the value of R can be lowered at will, down to 1/M, while keeping Tc fixed in order to 
generate the correct value for the four dimensional Plank scale from the 'low' fundamental 
scale ~ 1 TeV . To point out the importance of our result in this regime, we stress that our 
model makes possible the existence of a single extra dimension of very large size without 
contradiction with present observations, while keeping the value of the fundamental Plank 
scale around the TeV. 

The mechanism of suppression of heavy modes in our framework can relax the high 
energy constrains on the models with submillimeter extra dimensions. 

To summarize, there are two windows of the parameter space in the present model. The 
first one, in which the extra dimension has astronomic scale , is strongly constrained by the 
solar system observations which force R > 10^^ m (for R/rc — 10~^). In the second scenario 
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R can be anywhere below the Earth radius and 1/M. 



C. Star Cooling 

The most severe constraint on theories with low scale gravity comes from astrophysics 
since a stellar object can cool by radiating bulk gravitons 0. Let us estimate the rate of 
star cooling in the present case. Let the temperature inside the star be T. Then only the 
gravitons with masses m„ < T are produced efficiently. The rate of emission of the nth KK 
graviton is given by 

r„~r3G„. (3.16) 

We have to sum over all states up to n ~ TR. This gives 

2T^ R ^1 T^R aT^ 

This is less than the rate of just zero mode graviton production, and is totally negligible. 
Note that the contribution coming from the emission of a pair of KK gravitons via the virtual 
zero mode graviton exchange is of the same order, and is practically insensitive to the value 
of R and M. To see this let us consider a scattering process on the brane with a virtual 
zero mode graviton emission which subsequently decays into a pair of KK gravitons. The 
example of such a process can be a gravi-brehmstrahlung process, an electron scattering in 
the field of a nucleus and radiating gravitons. Let the typical energy in the process be E. 
For definiteness we take E > M, in order to put the strongest constraint. The rate of the 
process then becomes 

E^M^ E^ 
r.-^(Mfl).„^, (3.18) 

where we have cut-off graviton momenta by M. Surprisingly, this rate is similar to the 
one of single KK production, due to the suppression of direct KK production rate. Most 
importantly it has no dependence on either R or M. Thus the KK emission in high-energy 
processes a priory places no bound on either of this parameters. However, we do not expect 
M to be smaller than the inverse mm scale, due to no observed deviation from Newtonian 
gravity on larger distances, which gives the upper limit on the size of extra dimension 
R ~ alO^^cm. Thus the star cooling places no constraint on our scenario. It is useful to 
contrast this with the scenario of ^ . There, the dominant contribution to the cooling process 
comes from the production of heaviest modes, due to their high number. In the present case 
however, the wave functions of the heavy modes are suppressed on the brane and therefore 
cannot be produced efficiently. This is a generic property of the given framework with the 
induced kinetic term: in a high-energy process on the brane, the lightest KK states are 
the ones produced most efficiently. The extra dimension is more transparent for the softer 
modes! 
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D. Cosmology 



In this section we will consider cosmological constraints coming from the overproduction 
of bulk states. Not surprisingly, in analogy with star cooling, the cosmological constraints 
turn out to be rather mild. In order to be as model independent as possible, we shall discuss 
the following initial conditions for the hot big bang: 

1) The bulk is virtually empty; 

2) The brane states are in thermal equilibrium at some temperature Thrane, which can 
not be higher than the "normalcy" temperature T*. 

The normalcy temperature is defined as the temperature below which the Universe ex- 
pands as being effectively four- dimensional. This requirement then automatically restricts 
Tbrane < 7"*. Obviously, we would like X,, to come out at least as high as the nucleosynthesis 
temperature. As we shall see below, could be much higher, even as high as the cut-off scale 
(T* ~ M). This fact allows for the standard nucleosynthesis scenario to proceed unaffected, 
and also to implement one of the conventional baryogenesis mechanism, for generating the 
baryon asymmetry of the Universe. 

As discussed in [^, in theories with large extra dimensions the overproduction of bulk 
KK states can alter the standard cosmological expansion in two different ways. First, the 
energy density on the brane changes due to the "evaporation" into the bulk states. Secondly, 
the produced bulk states may dominate the energy density and over-close the universe. Let 
us consider constraints coming from these two effects separately. 

Cooling by Evaporation into Bulk States 

As estimated above for star cooling at temperature T the evaporation rate into the bulk 
gravitons is given by ( p.l7|) . The resulting change of the matter energy density on the brane 
due to evaporation is given by 



-^\evaporation ~ —T'^^ totaliT) ~ ' (3.19) 



Note that this is by a factor a smaller than the cooling rate due to the production of the 
standard zero-mode graviton, and is totally negligible in comparison to the cooling rate 
caused by the cosmological expansion 



^1 --xnn --,1! 

7 , expansion ~ "J-n p ~ n r 

at M, 



expansion 3H p ~ ~3-— — p , (3.20) 



where H is the Hubble parameter. For instance, in the radiation dominated epoch {H 
T'^/Mp), the ratio of the two rates is 



T 



dp I 

I evaporation 
{expansion j-vj.^ 



a. (3.21) 



This is a very small number even for T ~ TeV. The reason for this suppression can again 
be understood from the 'infrared transparency' of the theory: since the heavier gravitons' 
wave-functions are strongly suppressed on the brane, their production is not efficient enough 
to affect the brane-cooling process. 
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Overclosure by Gravitons 

Due to their suppressed couplings, KK gravitons are harder to produce on the brane. 
Moreover, due to the same suppression the KK modes are more long-hved as well. Thus, we 
have to examine constraints coming from their possible overproduction. 

The lifetime of a "cold" KK graviton of mass m can be written as follows: 

Tm . 3.22 

m 



This turns into the following relation upon using ( 1.19|) 



V M 



m 



(3.23) 



Thus, even the gravitons as heavy as m ~ M are stable for all the practical purposes and can 
over-close the universe if produced with a sufficiently large number density. This puts some 
bound on T* which we shall estimate below. Assuming that all the produced KK gravitons 
are stable and do not decay back to the brane, the energy density which is "pumped" into 
the bulk due to the brane-evaporation (in the process of the normal expansion) is 

Sp= -Trlevaporationdt ^ / ITJ-^C^r ~ IT^" • (3.24) 

Jti„ at JTeq Mp Mp 

Here T^g ~eV, is the standard crossover temperature and we neglected the later period of 
matter domination as well as the order one factors in the integral. Most of the energy density 
gets pumped in the bulk during one Hubble time after t = t^. Thus, soon after the initial 
time the bulk becomes populated by KK states with the energy density given by ( ^.241) . 



Most of these particles are relativistic and their energy density will redshift as radiation 
(p ~ T*^). However, to make our bound most conservative let us assume that most of the 
energy density is stored in the cold KK modes that redshift as matter (p ~ T^). Even in 
this case their energy density at the time of standard crossover will be given by 

PKK ~ ^aT^g . (3.25) 

Requiring that this be much smaller than the energy stored in the usual matter (at the same 
time t = teg), which, by definition, is pmatter ~ T^g ~ (lOeV)^, we get the following bound 
on T* 



Tl < ^eV = ^GeV' . (3.26) 



This bound can be easily satisfied even for T^, ~ M. 

In conclusion, we see that in contrast with the scenario of there are no essential cos- 
mological constraints from overproduction of bulk states due to the brane cooling. We have 
to stress, however, that this analysis can not capture more model dependent possibilities. 
For instance, if the KK states are produced by some other means in the early universe, 
they could either over-close the Universe or serve as unusual (and interesting) dark matter 



14 



candidates. To avoid the overproduction we have to assume that there was a period of the 
inflation that diluted the bulk and reheated only the brane 0. 



IV. ON THE ORIGIN OF THE CROSSOVER SCALE 

In this section we shall discuss a possible origin of the large distance scale Vc- The goal is 
to explain the large coefficient in front of the four- dimensional curvature term in Eq. ( p..8D . 
As already noted, in the effective field theory picture this term is not constrained by any 
symmetry and can emerge with an a priory unrestricted coefficient. However, it is desirable 
to have better understanding of this issue. As suggested in this term is induced due 
to quantum loops of the states which are localized on the brane and which interact with 
high-dimensional gravity. The resulting strength depends on the number of such states as 
well as their masses 0]. Thus, the large mass hierarchy could be obtained due to the large 
multiplicity of states propagating in matter loops on the brane. However, there can be other 
effects that may significantly contribute to the magnitude of this term. One possibility is to 
consider a brane-scalar field ^ which is non-minimally coupled to gravity 

^brane = / ^''x ^\detg\ CXp (^-^^ R{x) . (4.1) 

One could assume that, due to some dynamics, ^ develops a vacuum expectation value close 
to the scale M. This is in no contradiction with any fundamental principle, as long as ^ is 
somewhat lighter than M. Shifting the field —>■< ^ > +^ we end up with an exponentially 
large scale Vc = M~^e~^. Note that in the shifted vacuum the theory is consistent up to 
the energies of order ~ M. Moreover, the perturbative treatment is valid since the emission 
of ^ quanta are suppressed by the powers of M~^. 

Although at a first glance the fact that the VEV is larger than the cut-off of the theory 
might seem a bit unnatural, we should stress that there is nothing unusual in this fact. 
There exist many well defined examples both in string theory as well as in KK theories 
when this is the case. After all, the solution of the hierarchy problem in |jl|] can also be 
understood in this way. Indeed, the large value of the Planck scale is generated by the size 
of the extra dimension, or equivalently, by the VEV of the four-dimensional scalar field, the 
radion, which gets expectation value exceeding by many orders of magnitude the cut-off of 
the theory. 

Another well known example of this kind can be found in D-branes. It is well known that 
in the BPS limit the separation of two parallel D branes can be understood as the Higgs 
effect in the brane worldvolume gauge theory. The expectation value of the canonically 
normalized Higgs field is related to the string scale, Ms, as 

< e > ~ rM| , (4.2) 

where r is the inter-brane separation. This latter can be much larger than l/Ms- In fact 
the limit r » is well defined perturbatively and corresponds to the infrared limit 



^The inflationary solutions in the limit of a = oo were studied in |21] 
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for the bulk gravity. Although the world volume scalar field acquires an expectation value 
which is much larger than the fundamental scale, the theory is well defined. The two-brane 
system becomes populated with states of mass ~ rM| that are much heavier than the string 
scale (corresponding to the stretched string modes in the original theory). It may not be 
impossible that in a suitable framework like this, the brane separation can set the large 
coefficient of the brane-induced curvature term. 



V. IMPLICATIONS FOR SUB-MILLIMETER DIMENSIONS 

So far we have been studying the domain R < Vc- As we have shown, in such a case R 
can be of astronomical size. In fact, it can be arbitrarily large (if Vc is taken to be large). 
In this section, we shall study the opposite regime R > Vc- The case R ^ oo with one 
extra dimension was already discussed in 0. In this case the massive KK states which are 
lighter than Tc have unsuppressed wave-functions on the brane. Thus, they couple to the 
brane-matter with a strength comparable to that of the zero-mode graviton. Therefore, the 
masses of the KK modes must be large enough in order to avoid unacceptable deviations from 
Einstein's gravity. Precision measurements of light-bending and the precession of Mercury 
perihelion exclude the range of masses at the inverse astronomical scale. 

The shorter ranges are constrained by the precision measurements of G discussed above. 
These constraints push R to be somewhere around the inverse millimeter range, in fact. 



according to recent measurements |]I9| M~ > 250/i or so. In such a regime, planetary dy- 
namics is insensitive to extra dimensions, just like in the original scenario of However, the 
other predictions of this framework are dramatically modified. The situation is somewhat 
peculiar for <^ Vc ^ R. On the one hand, the predictions of []I| for table-top gravita- 



tional experiments |^ are unaffected. On the other hand all the current cosmological and 
astrophysical constraints are lifted, and the collider signatures are dramatically modified. 
This happens due to the suppression of heavy KK production at high energy colliders. 

In order to explicitly demonstrate this let us consider the star cooling process. This 
process puts the most stringent constraint on the scenario of for the sub-millimeter 
dimensions. For a given temperature T the production rate of nth KK state inside the star 
is suppressed as 

r~^V - (5 1) 



Assuming that rc ^ 1/T, we can evaluate the contribution of the KK modes with m < 1/r, 
and m > l/vc separately and show that these are both of order jgrf-- Using the relatioi 

Mp ~ M^i? we can bring the total contribution to the form 

rp4 



(5.2) 



Notice that there is an extra suppression factor ~ {Tr^ with respect to the standard case of 
This indicates that for sufficiently large Tc all the bounds can be avoided. Suppression 
of the bulk graviton production in other high energy processes can be analyzed in a similar 
way. The peculiarities of the spectrum of the model indicate that many of the experimental 



16 



constraints on theories with large extra dimensions must be reconsidered in the hght of the 
present discussion. 

VI. CONCLUSIONS 

In this paper we studied the scenario of Ref. with graviton kinetic term on the brane 
with one extra compact dimension. We showed that the existence of the graviton kinetic 
term on the brane allows for a novel framework with the high-dimensional fundamental 
Planck scale M ^ TeV (while the particle physics scale is bigger than a TeV), and a single 
flat extra dimension with the size which can be either smaller than 10^ m or larger than 10^^ 
m. 

The crucial role in generating usual 4D Einstein gravity on the brane is played by the 
brane-induced graviton kinetic term of reference 0]. The strength of this term is governed 
by a distance scale Tc which together with the size of the extra dimension defines the value 
of the four-dimensional Planck mass (see ( |1.19|) ). 

There are two phenomenologically interesting regimes. The first one is achieved when 
R Tc > 10^^ m. In this case the four-dimensional gravity on the brane is mediated 
by a single non-localized zero mode graviton, both at large and at short distances. The 
mass spacing of the KK modes is similar to that of an ordinary flat compact dimension. 
However, the wave-functions of heavy KK modes are suppressed on the brane by the ratio 
R/vc- This gives rise to the effect of "infrared transparency" |Q of the extra space. We 
studied the constraints imposed by precision gravitational measurements at all scales, as 
well as restrictions due to various astrophysical and cosmological effects. We found that the 
model is compatible with all those data. A crucial experimental test of this scenario could 
arise by observing deviations from Newtonian and Einstein's gravity practically at any scale. 
This includes precision studies of both relativistic and non-relativistic effects. In the present 
framework the states that live in the bulk can be a rather exotic source of dark matter. At 
distances r ^ R their interaction is that of an ordinary dark matter. However, for r <^ R 
they interact via much stronger gravitational potential which scales as Furthermore, 
these states interact with the observable matter by much weaker gravitational force. 

Another interesting limit is R> Tc and N > 2 [in D = 4 + N dimensional space). In this 
case the compactification radius R is constrained to be in a sub-millimeter domain. Although 
this seems to be similar to the scenario of [|l|, nevertheless, the framework is modified 
dramatically due to the graviton kinetic term on the brane. For instance, if < < R, 
the table-top predictions of scenario are unaffected and one should still expect deviations 
from Newton's law at scales r ~ i?. However, production of heavy KK gravitons is strongly 
suppressed. This lifts all the high-energy constraints. We have explicitly demonstrated this 
fact by reevaluating the constraints coming from the process of star cooling. Our analysis 
demonstrates the crucial importance of the brane-induced graviton kinetic term for the 
phenomenological studies. 

Note added 

After this work was completed, the related work ||22| appeared, which discusses the role 
of brane-induced kinetic term on compact dimensions in an example of a scalar field theory. 
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VII. APPENDIX 

In this Appendix we derive the expression for the force mediated by a scalar field in 
five-dimensional space-time with one compact dimension and an induced kinetic term on 
a 3-brane. We denote the coordinates x"^ = {x^,y),y G [0,-R]. Our starting point is the 
following five- dimensional Lagrangian 

C = M^d'^^dA^ + M'^S{y)d^'^d^^, (7.1) 

where ^{x^) has dimension 0, M is the fundamental five- dimensional mass scale, and 
is the scale which is dynamically generated by the interaction with matter that lives on the 
brane B . If we absorb the factor M'^/^ into the redefinition of the field then $ becomes a 
canonically normalized 5D scalar. The ratio Tc = M /M^ appears in the field equation as 
follows: 

d^dA^ + r^6{y)d''df,^ = 0. (7.2) 

Let us decompose the field $ into the following modes $ = J2m 'Pn{y)o'n{x^) ■ c"n(a;^) satisfy 
the four dimensional Klein-Gordon equation (9^5^ -|- mf^)an = 0, where m^'s are to be 
determined. The functions (pniv) set the profiles of the field in the fifth dimension. From 
( [7. 21) (with the use of d^dA = d^d^j, — dy) we get the "Schroedinger equation" for (pijj) 

[dl + m^ + r,mH{y)^(t){y) = 0. (7.3) 

This equation is to be considered on a circle of length R with the periodicity condition 
(f){y + R) = (piy). The problem is equivalent to that of a wave-equation for an infinite 
space with an array of delta-function type potentials located a.t y = nR, n = —0, ±1, ±2, . . ., 
with the identification y = y + R. The procedure for finding the solutions is standard: We 
solve the equation in two neighboring regions / = [— i?, 0], // = [0,i?], match the wave- 
function and its derivative at the boundary, and then identify the two regions by imposing 
the periodicity condition (f)j{y — R) = (f)jj{y). The solutions in each region are 

J) (f)j{y) = Ae'^^y + Be-'^^y 

11) (f)jj{y) = Ce'^^y + De-'™^ (7.4) 

where the coefficients A,B,C,D are to be determined by the following conditions: (a) 
periodicity; (b) continuity at y = 0; (c) matching of first derivatives at ?/ = 0: 
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b) A + 5 = C + D 

c) C - D - A{1 + zrcm) + 5(1 - irjn) = 0. 



(7.5) 



Solving these algebraic equations we determine the coefficients B, C, D in terms of A. This 
latter can in turn be found from the normalization condition on ch 



D = A B = C = A e"™^ , 

amended by the quantization condition for the masses 

Tc / ■mR\ / ■mR\ 

The nonlinear equation ( [7.7| ) is solved by graphical method as shown in Fig. 2 



(7.6) 



(7.7) 




Figure 2: Nonlinear condition for the mass spectrum (|7.7| ). The x axis is given in units of mR/2. 
The functions shown are tan and the hnes with slopes —Tc/R = —10, —1,0. For rc/R — the spectrum is 
standard Kaluza Klein spectrum m„ = 2mr/R, while for large rc/R it is given by m„ ~ (2n — 1)tt/R. 

We see that for R the masses approach the usual Kaluza-Klein spectrum, m„ = 

27m/_R. For R all modes (except zero mode) approach the asymptotes of the tangens: 
mo = 0, TUn — (2n — 1)71 /R. The level spacing for fixed R does not change. For arbitrary Tc 
the mass of the n'th state is in the interval ((2n — l)7r/R, 2mr/R). 

Let us now study the behavior of profiles in the fifth dimension. First we have to 
normalize (f)n{y) to unity. This fixes the coefficient An 



Ar 



An 



1 



rc/R 



l+rlml/4: 



2VR 



(7.8) 



For ^ R the coefficient A depends very weakly on the mass of the state and can be 
approximated by 1/ ^/R for all the modes. In Fig. 3 we show the modulus squared of 
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the wave-function for the zero mode (constant) and the lowest three modes for the choice 
rjR = 10. 



2.5 



1.5 



0.5 




Figure 3: Moduli squared of the wave functions for the zero mode and three lowest modes. The fifth 
coordinate is shown on the x axis and the units are R/2 (i.e. the whole range of the extra dimension is 
covered). The modulus squared is plotted on the y axis and the units are 1/R (i.e. the modulus of the zero 
mode is equal to one). The dependence is shown for the zero mode (constant), and lowest three modes for 
Vc/R — 10. The position of the brane at the origin is figuratively sketched as the thin "barrier" . 

The quantity which determines the couphng of the KK modes to the brane matter is the 
modulus squared of the wave-function at y = 0: 



i0n(o)r 



A + B\ 

2 



\A\ 



1 + e 



-im„R 



1 + tan^ (m„i?/2) 



\A\ 



1 + r1m\/A 



(7.9) 



Here, in the last line we used the condition ([7.7D. From ( |7.9| ) we deduce that higher KK 
modes are suppressed on the brane compared to the zero-mode. The suppression factor is 
exactly the same as in theories with infinite volume extra dimension |^. The mass of the 
n'th mode is of order 2?™//?, so we can rewrite ( |7.9| ) in terms of our initial parameters i?, Tc 



|0o(O)|' i + n^T^rjny 



(7.10) 



For Tc ^ R this suppression is substantial even for the lowest massive KK states. In 
Fig. 4 we show the modulus squared of the three lowest massive modes near the origin for 
rjR = 10. 
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-0.04 -0.02 0.02 0.04 

Figure 4: The same as Figure 3. The small portion of space near the origin is shown in order to 
emphasize the suppression of the first three massive modes on the brane. The zero mode has value one, 
which is well out of the y axis range shown in the figure. 

Let us now compute the potential due to the exchange of all the modes between two 



static sources located on the brane. For this purpose we add a source term to (fTTI), with a 
coupling to the field $ proportional to a/Gs = M~^/^. 

Decomposing the five dimensional field in four dimensional modes, $ = J2n4'n{y)o'n{x'^), 
and integrating (|7.1| ) over the extra compact dimension we obtain the following effective 
four dimensional Lagrangian: 

-^4 = 51 ^^CTm^/^CTn ( / dy (j)^(j)n ) + {(^m(^n) ( / dy 
n,m=0 \ V / V 

+rc (0™(O)0„(O)) + a„, . (7.11) 

Doing the integrals with respect to the compact coordinate dy one should take into account 
that the functions 0„ are not orthogonal. Nevertheless, we can proceed as follows: The field 
equation for 0^ can be written as: 

<Pm {dl + ml + r,ml6{y)) 0„ = 0. (7.12) 

On the other hand, 

0n {dy + r,ml6iy)) 0„ = 0. (7.13) 
Integrating both equations w.r.t. y and subtracting them, we obtain 

0={ml- ml) fdy{l + rj{y)) (7.14) 



This implies that the integral on the right hand side vanishes unless m = n. Using this fact, 
we obtain 

/ dy (t)ndy(t)m = ml dy(f)n4>^ (1 + S{y)rc) = ml6mn (l + rc|0„(O)|^) . (7.15) 
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Furthermore, inserting ( [7.15|) in (|7.11|) we get 



0(J 



n=0 



+ 



.(0) 



M3/2 



(7.16) 



In order to normalize canonically the kinetic terms we absorb the factor (1 + rc|(/)„(0)p)^/^ 
into the redefinition of the fields. Thus, the low energy Lagrangian becomes 



1 



0n(O) 



n=0 



y'l + r,|0„(O)|2VM^ 



CTnP 



(7.17) 



From (|7.17|) we can read off the couplings of the various modes to four dimensional matter. 
In the case of Tc ^ R, we can approximate m„ ~ (2n — 1)tt/R. Using this and some 
simplifying algebra we derive: 



Go 



Gn 



1 



Gr 



1 



R 



r^M^ \rj [{2n - l)7r/2]2 + R/r^ + {R/r^f 
R 1 



,M3 



(7.18) 



(7.19) 



An equivalent way to extract the information about the potential is to compute the 
propagator of the field $ from the Lagrangian ( |7.1| ). For this we have to solve the following 
equation (for simplicity we work in Euclidean space) 



^{x,y) = -5\x)5{y), 



(7.20) 



where A has dimension two. The procedure is similar to that used in with the only differ- 
ence being that the fifth dimension is now compact. We turn to 4-dimensional momentum 
space 



A(p,i/) = -Ky). 



(7.21) 



and take the ansatz A{p,y) = D{p,y)B{p), with D{p,y) satisfying 

{^dl-p')Dip,y) = -5iy). (7.22) 

This gives 

B(p) = -. (7.23) 

Eq. (|7.22|) can be solved in the compact space by expanding both sides in Fourier modes 



+00 -| +oa 

D{p,y)= E e^^^D.(p) S{y) = - ^ e^^^. 



(7.24) 
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As a result, one finds: 



1 



(7.25) 



To calculate the force between sources localized on the brane, we need to evaluate the 
propagator at ?/ = 0. Using the expression 



^coth^^^ 



2 / 



(7.26) 



we derive 



A(p,0) 



1 Tr. 



{pR/2) coth(pi?/2) 



p2 1 + {rc/R)ij)R/2) coth(pi?/2) 



1 



1 + (l/prc)2tanh(pi?/2) 



(7.27) 



The propagator never deviates substantially from over the whole range of p. The large-p 
and small-p behavior are respectively 



A(p,0) 
A(p,0) 



1 


^ 2 




pre 


1 


\-^] 


p2 



pi? ^ 1 (short distances) 
pR (large distances). 



(7.28) 



The maximal deviation in the coefficient of 1/p^ is of order R/re as before. If we continue 
this expression to Minkowskian space p —>■ ipM, the propagator becomes 



A(p,0) 



1 

pIi 



l + {l/pMrc)2tan{pMR/2) 



(7.29) 



The poles in this propagator are located at TcPm = — 2 tan(pM-R/2), in agreement with 
equation (|7.7| ) . 
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